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History 


Combustion  instability  is  an  organized, 
oscillatory  motion  in  a  combustion 
chamber  sustained  by  combustion. 

_ _ _ / 


\ 

Cl  caused  a  four  year  delay  in  the 
development  of  the  F-l  engine  used  in 
the  Apollo  program 

>  2000  full  scale  tests 

>  $400  million  for  propellants 
alone  (2010  prices) 

_ _ _ J 


Irreparable  damage  can  occur 
in  less  than  1  second. 


Damaged  engine  injector 
faceplate  caused  by  combustion 
instability 


"Combustion  instabilities  have  been  observed  in  almost  every  engine  development 
effort,  including  even  the  most  recent  development  programs" 

—  JANNAF  Stability  Panel  Draft  (2010) 
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Longitudinal  Experiment 


Continuously  Varying  Resonance  Chamber 
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Mean  pressure,  1.37  MPa  (200  psi) 
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Pulse  Timing  I 


Fuel  cut  off 


5.5"  Case 


Heat  release 
upstream  of 
the  backstep 


Flame  moves 
downstream 


•  In  the  3.5"  case  the  returning  pulse 
arrives  sooner  in  the  cycle 


In  the  7.5"  case  the  returning  pulse 
arrives  later  in  the  cycle 
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Pulse  Timing  II 


After  the  fuel  cutoff  event  the  combustion  restarts  through 

one  of  two  identified  mechanisms 

_ _ _ / 


Post-coupled  Ignition 
Mechanism:  The  returning 
wave  in  the  oxidizer  post 
pushes  unburnt  fuel  into  the 
warm  recirculating  gases  at  the 
backstep  where  ignition  takes 
place 


Vortex  Transport  Mechanism 

The  post  wave  arrives  later  in 
the  cycle,  the  unburnt  fuel 
slowly  mixes  with  the 
recirculating  gases  in  the  shear 
layer.  Ignition  takes  place  in  the 
shear  layer  downstream  of  the 
backstep 
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Baroclinic  Torque 
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MntheCVRC  ^ 
pressure  pulses 
in  the  oxidizer 
post  interact  with 
the  shear  layer 
and  generate 
\^vorticity _ 


Fuel 


Oxidizer  I 


Generation  of  vorticity  due  to  misaligned 
density  and  pressure  gradients 
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1 

Vp 
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Increased  vorticity  results  in  increased 
mixing  and  generates  partially  premixed 
regions  susceptible  to  combustion 
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Tribrachial  Flame 


•  Intersection  of: 

—  Diffusion  flame 

—  Partially  premixed  fuel 
rich  flame 

—  Partially  premixed  fuel 
lean  flame 

•  Location  of  intense  heat 
release 

•  Observed  in  the  CVRC  by  Garby 
et.  al.  and  Guezennec  et.  al 


White  line,  T  =  2000  K 
Black  line,  Z  =  Zsl  =  0.095739 
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2D  Parametric  Study 


2D  Swirl  shows  substantial 
reductions  in  amplitude 
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Concerns: 

Total  amplitudes  in  2D  are  lower 
compared  with  3D 


-  Swirl  component  is  assumed 
axisymmetric  and  may  be 
stronger  than  3D 
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2D  Shows  Attached  Flame 
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Three-dimensional  Study 


Examine  the  effect  of  swirl  when  exposed  to  the  larger  amplitude 
consistent  with  the  experimental  results 


Remove  the  axisymmetric  assumption 


Swirl  is  imposed  using  a  boundary  condition,  not  geometry 


uq  (r)  =  u 
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R. 


sin(6»s) 


Four  cases: 


Baseline 


3°  Swirl 


9°  Swirl 


15°  Swirl 
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Simulation  without  Swirl 


Mach  Number 


25%  CH4  Temperature,  K 

400  800  1200  1600  2000  2400  2800 

# 

Heat  Release  Surface 
Colored  With  Pressure 


1.4E+06  1.5E+06  1.6E+06  1.7E+06  1.8E+06  1.9E+06 
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Simulation  with  Swirl 
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Pressure  Amplitude 


(  \ 

Swirl  has  a  wider  trough  with 

more  unsteady  fluctuations 
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PSD  Power,  kPa2/Hz  PSD  Power,  kPa2/Hz 


Power  Spectral  Density  Analysis 
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Amplitudes 


\ 

PSD  data  is  integrated  using  the  FWHM  method  to  determine  the  peak-to- 

peak  pressure  for  each  mode. 

_ _ _ _ _ J 


Mode 
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/,  Hz  p\  kPa 

/,  Hz 

3° 
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/,  Hz 

9° 
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2 
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3 
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9°  of  swirl  had  the  greatest  reduction,  about  30% 


Amplitude  reduction  is  less  than  the  2D  prediction 


The  unknown  mode  has  a  weak  amplitude. 
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Methane  Cycle 


(a)  Time  1. 


CHq  Mass  Fraction 
Static  Pressure,  MPa 


0.05  0.15  0  35  0,35  0.45  0,55  0  85  0  75  0  35  0  95 


1.20  1 .26  1.32  1.38  1.44  1.50  1.56  1  62  T  .68  1.74  1.80 


(c)  Time  2. 
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(d)  Time  2. 


(f)  Time  3. 


9°  Swirl 


17 


Methane  Cycle  II 


(g)  Time  4. 


(h)  Time  4. 


(i)  Time  5. 


(j)  Time  5. 
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Heat  Release  Cycle 


(a)  Time  1. 


(b)  Time  1. 


(c)  Time  2. 
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Heat  Release  Cycle  II 


(g)  Time  4. 


(h)  Time  4. 
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Time  6. 


9°  Swirl 
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Lack  of  a  CTRZ 


Central  toroidal  recirculation  zones  can  be  found  in  flows  with  swirl. 

\ _ J 

The  CTRZ  has  a  stabilizing  effect 

V _ ) 


f  \ 

A  CTRZ  was  found  in  2D  but  not  3D 

V _ ) 


; 

Possibly  a  result  of 
the  artificial 
centerline  boundary 
condition 

_ _ _ ) 


Approved  for  public  release;  distribution  unlimited 


21 


Time  Averaged  Flowfield 


Heat  Release, 
MW/m3 


2000  8000  14000  20000  26000  32000  38000  44000  50000 


Heat  Release, 
MW/m3 


2000  8000  14000  20000  26000  32000  38000  44000  50000 


Heat  release  is  similar 
for  all  three  swirl  cases 
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K. 


Swirl  shows  that  the 
heat  release  that  is 
confined  to  the  shear 
layer 
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J 


No  "bump"  in  heat 
release 


No  evidence  of  a  CTRZ 
in  3D! 
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Accumulation  of  Fuel 


Cycle  analysis  showed  that  heat  release  resumes  quickly  in  the  swirl  case 

compared  to  the  baseline 


v_ 


There  is  less  fuel  accumulation  in  the  cycle  and  evidence  of  consumption 

sooner 
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Injection  Differences 


Temperature,  K  Temperature,  K 

300  560  820  1080  1340  1600  1860  2120  2380  2640  2900  300  560  820  1080  1340  1600  1860  2120  2380  2640  2900 


Baseline 
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Injection  Differences 


Baseline 


More  radial 
expanse 


Less  pure  fuel 
regions 
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Injection  Differences 


More  radial 
expanse 
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Injection  Differences 


More  radial 
expanse 


Less  evidence  of 
fuel 
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Pressure  Amplitude,  Pa 


Unknown  Mode 
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Unknown  mode  was  found  in  the  swirl  cases  with  a  frequency 

between  the  first  and  second  longitudinal  modes 
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(a)  3°  swirl,  pressure. 


(b)  9°  swirl,  pressure. 


(c)  15°  swirl,  pressure. 


\ 

No  evidence  of  a  PVC  mode  or  spinning  mode 

_ J 


"  ^ 

A  similar  mode  was  found  in  a  3.5"  simulation  which  had  a 

similar  amplitude  to  the  swirl  cases. 

_ _ _ J 
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Summary 


A  numerical  investigation  of  the  effect  of  swirl  on  an  unstable  single-element 

gas-gas  rocket  engine  was  undertaken 


The  addition  of  swirl  to  the  fuel  reduced  the  amplitude  by  30% 


The  reduction  of  amplitude  was  a  result  of: 


Improved  mixing 
between  the  fuel 
and  oxidizer 


Shorter  recovery 
time  of  the  fuel 
cut  off  event 


More  continuous 
heat  release 
throughout  the 
cycle 


The  percent  amplitude  reduction  in  3D  was  not  as  great  as  2D,  but  had  a  similar 

absolute  reduction  in  amplitude 


2D  simulations  should  be  used  to  reduce  the  number  3D  simulations 
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